Mass and Cross-section Measurements of 
Chargino at Linear Colliders in Large tan (3 

Case 



Yukihiro Kato a '*, Keisuke Fujii b , Teruki Kamon , 
Vadim Khotilovich c , Mihoko M. Nojiri d 

^Department of Physics, Kinki University, Osaka, Japan 
h High Energy Accelerator Research Organization, Ibaraki, Japan 
c Department of Physics, Texas A &M University, College Station, TX, U.S.A. 
d Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto, Japan 



Abstract 

The lighter chargino pair production (e + e~ — ► XiXi) 1S one °f t ne key processes 
for determination of supersymmetric parameters at a linear collider. If tan/3 is a 
large value, the lighter stau (t^) might be lighter than the Xii while the other 
sleptons stay heavier. This case leads to a cascade decay, X\ ~* t~\ v t followed by 
— > Xi r± ) where Xi is the lightest supersymmetric particle. This paper addresses 
to what extent this cascade decay may affect the measurements of the chargino 
mass and its production cross section. 
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1 Introduction 



Supersymmetry (SUSY) is one of the most promising scenarios of physics be- 
yond the Standard Model (SM) and predicts new particles observable at the 
present or next generation of colliders. It solves the naturalness problem by 
taming the otherwise quadratically divergent quantum correction to the Higgs 
boson mass, makes the three gauge forces unify, and opens up a way to ultimate 
unification of all the fundamental interactions including gravity |l| . If SUSY is 
the case, we will be able to probe physics at the unification scale by studying 
its breaking pattern through various measurements of SUSY particles some 
of which will have been discovered at the Tevatron or the LHC One of pri- 
mary goals for an International Linear Collider (ILC) is hence determination 
of SUSY or SUSY-breaking parameters through precision measurements of 
masses and couplings of those SUSY particles. The pair production of lighter 
charginos (e + e~ — > xtxi) is one °f the most important processes in this re- 
spect since the xt is relatively light in most of the parameter space and hence 
expected to be studied in the early stage of the ILC project. Measurements of 
its mass and couplings at the ILC have not only been discussed but, for some 
typical sets of SUSY (-breaking) parameters, also simulated in detail 0. 

All of the cases studied in the past, however, assume that the lighter chargino 
would decay directly into a final state containing the lightest supersymmetric 
particle (LSP), which is the lightest neutralino (xt) in most of the SUSY 
parameter space. At large tan (3 values, however, tl-tr mixing might generate 
a large mass splitting between f\ and the other sleptons (£ = e, jl, f 2 , u), 
yielding a particular mass hierarchy of M~ t > M~± > M~± > M^o. If xt 
is gaugino-like, this means that the xt decay at 100% of the times in 
the cascade mode: xt ~^ ^i u t followed by ff 1 — > Xi T± - Consequently, the 
experimental signature of the chargino pair production will be an acoplanar r 
pair with significant missing energy carried away not only by the LSP's from 
the fi decays but also by the final-state neutrinos from the xt an d t decays. 
This would make it significantly more difficult to determine the chargino mass. 

This paper describes results from our simulation studies of the pair production 
of the lighter charginos and their subsequent cascade decays: e + e~ — > XiXi ~^ 
f 1 f z/ r + f 1 ~z7 T in a particular scenario of the Minimal Supersymmetric Standard 
Model (MSSM), where M-± = 152.7 GeV/c 2 and M^o = 86.4 GeV/c 2 as in 
Refs. Hfl. It should be noted that signals from the f^xt co-annihilation 
region are recently of great interest. The region is characterized by a differ- 
ence in the mass and the lightest neutralino mass of about 5-15 GeV/c 2 . 
This small mass difference allows the f ± to co-annihilate in the early universe 
along with the neutralinos in order to produce the current amount of cold 
dark matter density of the universe measured by WMAP 0. However, the 
experimental signature is very challenging. We will not address this signature, 
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because it requires to detect electron/positron from two-photon process in the 
forward region. Thus it requires a different optimization of the event selection 
criteria 0,13, - 



2 Monte Carlo Simulation 



In simulating the signal and various background events, we use the HELAS 
library |9J to calculate the helicity amplitudes, the BASES /SPRING pack- 



age jlOj] to generate the final-state partons, and TAUOLA to decay the 
final-state r leptons if any, so as to take into account possible spin correla- 
tions of intermediate heavy partons, effects of beamstrahlung and subsequent 
beam-energy spread and finite r polarizations. The generated events are 
then processed through a fast simulator of a detector model 0], which in- 
cludes a central tracker, electromagnetic and hadron calorimeters, and muon 
drift chambers. This simulator smears charged-track parameters in the central 
tracker with parameter correlation properly taken into account, and simulates 
calorimeter signals as from individual segments. Although r polarization varies 
in general with the SUSY parameters and would thus influence the property of 
the final states, we fix the r polarization for the signal process at the nominal 
value of 0.026 for simplicity. 

For SUSY events, we choose our reference SUSY point to have tan/3 = 50, 
[i = 400 GeV, m = 200 GeV, and M 2 = 180.0 GeV, where //, m , and 
M2 are the higgsino mixing mass, universal scalar mass at GUT scale, and 
SU(2) L gaugino mass parameters, respectively. We take Mi so that it satisfies 
GUT relation M 1 = 0.5M 2 . It corresponds to M~± = 171.0 GeV/c 2 (M r j£), 

M^o = 86.4 GeV/c 2 , and M-± = 152.7 GeV/c 2 . Notice that ff is the next to 
the lightest SUSY particle (NLSP) and decays to jfir at 100% of the time. We 
assume in what follows that M^o and M-± will have been measured below the 
threshold of chargino pair production, say at a/s = 310 GeV, with uncertainties 
expected for / Cdt = 100 fb _1 |4J- In order to evaluate to what extent we 
can determine the chargino mass in this setting, we then vary the Xi mass 
around 171.0 GeV/c 2 , while keeping the ff* 1 and Xi masses unchanged, by 
appropriately adjusting the values of SUSY(-breaking) parameters 1 and see 
how the observable distributions respond to the change. 

As for the ILC parameters, we assume y/s = 400 GeV, / Cdt = 200 fb _1 , 
and an electron beam polarization of V(e~) = —0.9 (left handed). The X1X1 
production cross section at V(e~) = —0.9 is 20 times larger than that at 



1 For this we allow violation of the GUT relation between the U(l) and SU(2)l 
gaugino mass parameters (Mi and M2) that has been assumed for our reference 
point. 
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V(e~) = +0.9 (right handed). Furthermore, it can suppress the cross section 
for the e + e~ — > f{~ process that yields almost the same final states as our 
signal process. 



3 Event Selection 



In our setting, the chargino pair production results in a final state consisting 
of a r-pair with large missing energy. The five major backgrounds for the sig- 
nal event are considered: (i) f^-pair production: f/ff — > t + x\ + t — Xi? ( n ) 
neutralino pair production: X1X2 ~^ Xi + ~~ > + r + r~, (iii) diboson 
production: W + W~ -> r+z/ r + r"F T , -> t + t~ + (£ = e,/i,r) or 

r + r~ + z7i/, (iv) single boson production: eeZ° and vvZ° with Z° — > £ + £~ 
or z/77, and (v) two-photon r-pair production: e + e~ — ► e + e~r + T _ . As for the 
neutralino pair production, the SU(2)l relation between X2 an d Xi masses is 
assumed. We also assume the X2 m ass is measured via other processes such 
as X2X2 ~* 4t + 2Xi- We ignore the uncertainty in the X2 m ass measurement, 
because the effect in Xi mass measurement is found to be negligible. 2 Notice 
that the cross section times branching fraction for 1^-pair production followed 
by W — > iv decays is, for instance, eight times higher than that of the signal 
at our SUSY reference point. We had thus better avoid events with purely lep- 
tonic r decays. Notice also that the two r-jets 3 in the SUSY events are likely 
to be more acoplanar than those from the SM backgrounds. These observa- 
tions lead us to the following selection criteria to improve signal-to-background 
ratio: 



(a) there should be no e/fi candidates positively identified in the detector 
simulator; 

(b) jet clustering has to yield two jets with 5 GeV < £j et < 160 GeV for 
Vent > 0.0025 Q; 

(c) — Qi ■ Ji Z /\Ji\ < 0.8, where Qi and Jj are the charge and the momentum 
of i-th jet with % — 1(2) corresponding to the higher (lower) energy jet; 

(d) M jet < 3 GeV/c 2 ; 

(e) the missing transverse momentum P™ lss > 20 GeV/c; 

(f) cos0(Ji, P v i S ) < 0.9 or cos9(J 2 , P v is) > —0.7, where cos 6( Ji, P v i S ) = Ji ■ 
-Fvis/d^ill-FvisI) an d -P v is is the visible momentum vector calculated with 
both tracker and calorimeter information; 

(g) thrust < 0.98; 

(h) acoplanarity angle 9 a > 30°. 



2 The x^ mass determination accuracy is estimated to be ±0.02 GeV/c 2 if the x\ 
mass fluctuation is ±1 GeV/c 2 . 

3 A group of particles coming from a r decay will hereafter be called a r-jet even 
if it contains only one charged particle. 
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Cut (f) effectively rejects the WW background, because the r lepton from 
W — > tv tends to go in the parent W-boson momentum direction. 

We tabulate the resultant selection efficiencies after each cut for the signal 
and the background processes at the reference SUSY point in Table 1. Total 
selection efficiency (e tot ) for the XiXi signal at the reference point is e[^ t = 
(15.8 ± 0.1)%. We see that the selection criteria are effective at suppressing 
the SM backgrounds, yielding a signal-to-background ratio of 1.7. As for the 
ZZ — > vvqq event where quark jets are mis-identified as r-jets, the total 
selection efficiency is small (less than 5.0 x 10 -5 ) and the expected number 
of the ZZ — > vvqq events is 2.5. We do not include this process in Table 1. 
We also show the signal selection efficiencies for different chargino masses in 
Table 2. The and Xi masses are fixed to be 152.7 GeV/c 2 and 86.4 GeV/c 2 , 
respectively. The statistical uncertainty in the efficiency for each point is better 
than 0.1%. The efficiency depends on the Xi mass only very weakly as long 
as AM = M~± — M~± is large enough to deposit a jet energy of E-, ct > 5 GeV 
m the detector. 



4 Results 

4-1 Determination of chargino mass 

If Xi an d T"i masses are predetermined, the r-jet energy distribution can be 
used to determine the xf mass, since the shape of the distribution, especially 
around the endpoint, depends on AM. As AM increases, the peak of the 
distribution shifts to the higher energy side, while the endpoint moves to the 
lower energy side. This can be seen in Figure 1, which plots the r-jet energy 
distributions for various AM values. 

In order to determine the chargino mass using the r-jet energy distribution, 
we first prepare a parametrized r-jet energy distribution (a template) as a 
function of M~± in the following manner: (i) the r-jet energy distribution, 
including the backgrounds, is fit to a polynomial function by MINUIT ^5| 
for each of high statistics samples with different Xi masses; (ii) the fitting 
parameters are then parametrized as a function of M^±. Although the r-jet 
energy distribution may depend on the remaining free parameters such as z> 
mass, \i parameter, and r polarization, those are fixed in the present study. 

We then compare the r-jet energy distribution for a Monte Carlo sample 
corresponding to 200 fb _1 with the template and calculate x 2 ■ Figure 2 shows 
Ax 2 = X 2 - Xmin for the reference point (Mi e £ = 171.0 GeV/c 2 ). The fitted 

Xi 

value of the xf mass is 170.5 ± 0.6 GeV/c 2 . This is about la off the reference 
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Table 1. Cumulative selection efficiencies (e) for signal (xtxi ) an d various sources of backgrounds. The cross section includes branching 
ratios of W — > eu,^u,Tu for W + W~ process and Z° — > ee, ^fi,TT,uu for Z°Z°, eeZ°, and z/zvZ processes. 



process 


(ref.) 


(ref.) 


(ref.) 


w+w~ 


Z°Z° 


eeZ° 


isisZ 


e + e — > 

e + e~T + r~ 


a (fb) (P(e-) = -0.9) 


322 


78 


44 


2304 


83.6 


1876 


92.8 


5.62xl0 4 


(a) e//x removal 


0.403 


0.424 


0.398 


0.0747 


0.0623 


3.7 x 10~ 2 


0.0437 


7.3 x 10- 1 


(b) N jet = 2 


0.351 


0.308 


0.348 


0.0392 


0.0522 


1.8 x 10~ 2 


0.0377 


2.6 x 10" 2 


(c) Jet angle cut 


0.240 


0.218 


0.251 


0.0122 


0.0255 


6.9 x 10~ 3 


0.0239 


2.4 x 10~ 2 


(d) M jct < 3 GeV/c 2 


0.220 


0.187 


0.238 


0.0118 


0.0229 


5.8 x 10~ 3 


0.0232 


2.4 x 10~ 2 


(e) P^ iss > 20 GeV/c 


0.182 


0.147 


0.198 


0.0085 


0.0208 


4.8 x 10" 4 


0.0178 


1.9 x 10" 4 


(f) cos0(Ji,P vis ) < 0.9 or 


















cos#(J 2 ,P vis ) > -0.7 


0.175 


0.142 


0.185 


0.0054 


0.0200 


4.7 x 10" 4 


0.0161 


1.2 x 10" 4 


(g) Thrust < 0.98 


0.173 


0.141 


0.172 


0.0043 


0.0200 


4.1 x br 4 


0.0160 


1.1 x 10~ 4 


(h) Acoplanarity 9 A > 30° 


0.158 


0.137 


0.158 


0.0031 


0.0186 


9.1 x HT 5 


0.0136 


2.8 x 10" 5 


Ctot 


0.158 


0.137 


0.158 


0.0031 


0.0186 


9.1 x 10~ 5 


0.0136 


2.8 x 10~ 5 


(generated events) 


0.4M 


0.4M 


0.6M 


1M 


1M 


1.4M 


1M 


8M 


# of expected events 


















for = 200 fb" 1 


10175 


2137 


1390 


1428 


311 


34 


252 


315 



M-± = 171.0 GeV/c 2 , M f ± = 152.7 GeV/c 2 



Table 2 

Total selection efficiencies for five different chargino masses where M~± = 152.7 
GeV/c 2 and M^o = 86.4 GeV/c 2 . 



M-± (GeV/c 2 ) 


161.0 


166.0 


171.0 
(ref.) 


176.0 


181.0 


AM(= M-± - M~±) (GeV/c 2 ) 


8.3 


13.3 


18.3 


23.3 


28.3 


a(fb) 


433 


378 


322 


268 


214 


^tot 


0.156 


0.157 


0.158 


0.158 


0.158 



pt = 400 GeV, tan/3 = 50, r polarization = 0.026, V(e~) = -0.9 




T-jet energy (GeV) 

Fig. 1. Energy distributions of r-jets in the signal + backgrounds events for 3 
representative choices of AM after our final selection criteria (Table 1). We also 
show background events. The shape in the reference SUSY scenario (AM = 18.3 
GeV/c 2 ) is compared to other cases (AM = 8.3 and 28.3 GeV/c 2 ). 

mass of 171.0 GeV/c 2 . We prepare approximately 100 statistically independent 
samples of XtXt an d X2X1 events for the reference point. Each sample is also 
corresponding to 200 fb -1 . The x 2 fitting procedure provides consistently a la 
uncertainty of ±0.6 GeV/c 2 for each sample. However the fluctuation of xt 
mass at the Xmin i s ±0.4 GeV/c 2 around the Ml±. Therefore, we assign the 

Xi 

total uncertainty in the xt mass measurement to be ±0.7 GeV/c 2 (quadratic 
sum of two sources). 

The above estimate ignores the uncertainties in the input fi and Xi masses. 
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liii.iiii.ii.n.i.*.«i.^.i.i.ii..i.ii.i.i. 
168.5 169 169.5 170 170.5 171 171.5 172 172.5 

Mj^(GeV/c 2 ) 



Fig. 2. Ax 2 of r-jet energy distributions for the reference point 

These uncertainties will affect the accuracy in the Xi mass measurement, 
because the differences between f x , Xii an d Xi masses determine the shape of 
r-jet energy distribution. The higher endpoint, in particular, depends on the 
difference (AM) between ff 2 and Xi masses. In order to evaluate the impact 
of the uncertainties on the xf mass determination, we repeat the study for 
M-± = 171.0 GeV/c 2 , varying M~± and M^o. Figure 3 shows a contour plot of 
Ax 2 in the M-± - M^o plane. The center point of each A% 2 contour is at M^o = 
84.6 GeV/c 2 and M~± = 151.2 GeV/c 2 . The contours clearly demonstrate a 
strong positive correlation between M^o and M~± , implying that the shape of 
the r-jet energy distribution is mainly controlled by their difference, as long 
as the xf mass is fixed. It is worth noting that the smallest ellipse, which 
is the expected la bound from the rf-pair production study with 100 fb _1 
cited from Ref. has a quite similar correlation and is well contained in the 
contour corresponding to Ax 2 = 1. The aforementioned expected accuracy 
of ±0.7 GeV/c 2 for the Xi mass will hence essentially be unaffected by the 
inclusion of the uncertainties in the Xi an d ri masses from the lower energy 
fj^-pair study. 

4-2 Production cross section measurement 

The production cross section times branching ratio is experimentally estimated 
by a ■ Br = (N obs — N bg )/(e tot ■ / Cdt), where N obs is the number of observed 
events, N bg the number of expected background events, / Cdt the integrated 
luminosity, and e to t the selection efficiency. It should be noted that e to t weakly 
depends on M~± as shown in Table 2: the efficiency varies by only 0.6% over 
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Fig. 3. Contour plot in the M~± — M^o plane, where M~± can be measured 

to an accuracy of ±0.7 GeV/c 2 (A X 2 = 1), ±1.4 GeV/c 2 (A X 2 = 4), and 
±2.1 GeV/c 2 (Ax 2 = 9). The smallest ellipse is the expected la bound for the 
Xi and mass fit with 100 fb -1 cited from Ref. 0|. 

the mass range of ±20 GeV/c 2 from Ml±. This implies the possibility of 

Xi 

extracting the cross section times branching fraction (er ■ Br = 322 ± 6 fb) 
independently of the input xt mass, provided that / Cdt = 200 ± 2 fb -1 
(the luminosity is assumed to be measured with an uncertainty of 1% 0]), 
Nobs — Nb g = 10175 events (the background contribution is fully estimated), 
and e to t = 0.158 ±0.002 (the uncertainty includes the weak dependence of e to t 
on M~±). This indicates that we can determine the production cross section 
times branching fraction to an accuracy of 2% even without any precise M^± 
study. 



5 Summary 



At large tan/3 values, the f{ mass might be lighter than the Xi rnass within 
the MSSM framework. In this case, the xf — > t^v t decay will be dominant 
and that will be followed by rf — > r ± x'i- We have thus studied the final state 
consisting of rr ± large missing energy as a possible signature of the chargino 
pair production at a future e + e~ linear collider. 

We proposed to use the r-jet energy distribution to measure the xt mass in the 
large tan/3 scenario. Since this scenario leads xi an d ff 2 masses to be lighter 
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than Xi mass and opens the f x pair production mode with lower energy of Xi 
pair production, the Xi an d masses are likely to be determined in the lower 
xf production threshold energy operation of the collider. The measurements 
of 171 GeV/c 2 chargino could then be made with an accuracy better than 
1 GeV/c 2 , provided that the Xi an d masses will have been determined 
through a f^-pair production study with 100 fb _1 . 

Since the event selection efficiency does not significantly depend on the Xi 
mass, the uncertainty in the production cross section shall be dominated by 
uncertainties in the selection efficiency and the luminosity measurement. We 
expect to determine the cross section to an accuracy of 2% at the reference 
point despite the cascade decay mode. 
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